A variational description of the quantum phase transition in the sub-Ohmic spin-boson model 
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The sub-ohmic spin-boson model is known to possess a novel quantum phase transition at zero temperature 
between a localised and delocalised phase. We present here an analytical theory based on a variational ansatz 
for the ground state, which describes a continuous localization transition with mean-field exponents for < s < 
0.5. Our results for the critical properties show good quantitiative agreement with previous numerical results, 
and we present a detailed description of all the spin observables as the system passes through the transition. 
Analysing the ansatz itself, we give an intuitive microscopic description of the transition in terms of the changing 
congelations between the system and bath, and show that it is always accompanied by a divergence of the low- 
frequency boson occupations. The possible relevance of this divergence for some numerical approaches to this 
problem is discussed and illustrated by looking at the ground state obtained using density matrix renormalisation 
group methods. 
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The physics of quantum systems in contact with environ- 
mental degrees of freedom plays a fundamental role in many 
areas of physics, chemsitry and biology, including systems 
as diverse as solid state quantum computerslJl El, quantum 
impurities [3], and photosynthetic biomolecules||4|-l8|]. A key 
theoretical model for the study of system-environment inter- 
actions is the spin-boson model (SBM), which consists of a 
two-level system (TLS) that is linearly coupled to an 'environ- 
ment' of harmonic oscillators JsKl^l- Although this model has 
been studied extensively, there are still many open problems 
in SBM physics, most notably those concerning the quan- 
tum phase transition (QPT) between delocalised and localised 
phases that exists in the SBM when the oscillators are charac- 
terised by sub-Ohmic spectral densities. 

The standard quantum-classical mapping predicts that the 
sub-ohmic SBM should be equivalent to a classical Ising spin 
chain with long-range interactions, and predicts a continu- 
ous magnetic transition with mean-field critical exponents for 
< s < 0.5. In Ref. lllll . a continuous transition in the 
sub-Ohmic SBM was observed using the numerical renormal- 
isation group (NRG) technique for all values of < s < 1, 
but the critical properties of the transition were found to be 
non-mean-field for < s < 0.5. It was suggested that this im- 
plied a breakdown of the classical to quantum mapping, and 
some subsequent work in this and other systems has supported 
this claim lll2Ul3[l . However, it is now believed that the non- 
mean-field results found by NRG in < s < 0.5 are incor- 
rect, and arise from the truncation of the number of states Ni, 
used to describe each oscillator in the Wilson chain lfl4l[l5ll . 
Recent studies of the sub-Ohmic QPT using quantum monte 
carlo (QMC) lll6ll . sparse polynomial space approach (SPSA) 
1 1711 . and an extended coherent state technique have indeed 
found mean-field critical exponents for < s < 0.5 lIlSll . 

In this article we propose a variational ansatz for the ground 
state of the sub-Ohmic SBM for < s < 0.5 which does 
not require any truncation of the environment. This is an im- 
portant feature, as we shall show that the number of environ- 
mental bosons diverges above the transition. Although our 
ansatz is essentially just a guess at the form of the ground 
state, we shall show that our results agree extremely well with 



existing numerical results and support the mean-field picture 
of the transition for < s < 0.5. Moreover, we are able 
to give an inuitive picture of the phase transition that might 
help in improving numerical simulations of the ground state, 
and illustrate this with a density matrix renormalisation group 
simulation of the SBM ground state. 



THE VARIATIONAL ANSATZ 
The spin-boson Hamiltonian can be written (h ~ 1) as 



H = Ada; 



^(Tz'^giiai + aj) + '^ujia'lai. (1) 



where <Ti are the usual Pauli operators which describe the TLS 
and ai , aj are the bosonic anihilation and creation operators 
respectively of bath modes of frequency oji. The tunneling 
amplitude of the TLS is A, and gi are the couplings between 
the TLS and the bath modes. It is well established that all ef- 
fects of the bath on the reduced state of the TLS are completed 
determined by the spectral function J(w) = tt J^i di^i^^^i) 
jgi fioll . Following Bulla et al. we consider the spectral func- 
tion J{lj) = 2Tr aijjl'''uj^O{ujc — i^) ll20ll . where uJc is the 
maximum frequency in the bath. Super-Ohmic baths have 
s > 1, Ohmic baths s = 1 and sub-Ohmic baths s < 1. 

Representing the +1, —1 eigenstates of az as |+), |— ) re- 
spectively, a variational ground state ansatz j^*) is written in 
the following way, 



1*) (2) 
where \4>±) is given by 

|0±> = e-^'^'±('^'-'^?)|0), (3) 

and |0) is the vacuum of the bath modes. This ansatz de- 
scribes a dressed TLS, a superposition of the localised states 
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|±) which are correlated with bath modes displaced by fi±. 
For general displacements, the wave function is not separa- 
ble into a product state of the TLS and environment, and 
our ansatz is not simply a mean field theory like the that re- 
cently studied by Hou et al. jlSll . Our ansatz is a generali- 
sation of the variational wavefunction of Silbey and Harris 
(SH) iST], in which the constants and displacements are fixed 
to obey C+ — C_ and /;+ = —fi-- As we shall show, 
these constraints are broken in the localised phase, and we 
shall henceforth refer to the ansatz of Eq. (|2) as the Asym- 
metrically Displaced Oscillator (ADO) state. The normali- 
sation of the wavefunction (^'|^') = 1 enforces the relation 
+ = 1, and we assume that C+, C_ are real. The or- 
der parameter of the TLS localisation transition will be taken 
as the magnetisation M = (^'Icr^ j^*), which can be expressed 
as M = — . Using the normalisation condition we ob- 
tain C|. = i(l ± -A/). With these relations and the standard 
properties of displaced oscillators, one obtains the following 
expression for the ground state energy E{M) = 



■^{fi+9i + fi+^i) 



(1 



] Z^Ui-9i - Ji 

I 



(4) 



In Eq. (IDi we have introduced a renormahzed tunneling 
amplitude A given by 



A 



Aexp 



(5) 



The quantity A/ A is the overlap of the displaced oscillator 
wave functions which dress the TLS states |±). We now min- 
imize the energy at constant ]\I w.rt. the displacements to 
obtain. 



.fi± = - 



gi{MA±Vl-M^uJi) 



(6) 



2w;(A + Vl -APuji) 

These displacements are then substituted back into Eq.©, and 
the sums converted into integrals using the spectral function. 
The integrals can be computed exactly in terms of hyperge- 
ometric functions, but in the scaling limit ujc co, these 
hypergeometric functions can be expanded in powers of w^^. 
For s < 0.5, the ground state energy to leading order in — 
takes the simple form. 



2s 

anujcil - s){l - M2) 
2 sin(7rs) 



A 



,V1 - 



(7) 



The renormalised tunneling amplitude is the solution of the 
implicit equation. 



Aexp 



(1 - AP)io''duj 
(ATTT^^^ 



(8) 




E(M)-E(0) 




FIG. 1. Ground state energy difference E{M) — E{0) as a function 
of magnetization for s — 0.3, A = 1 and coc ~ 10. The curves 
correspond to a = 0.325, 0.3265, 0.327 and 0.328 (black,green, red 
and blue lines respectively). 
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FIG. 2. Analytical critical coupling ac as a function of s for A = 1 
and various oJc. 



which in the scaling Umit takes the form 

A = A exp 



l-s 



1 — s sin(7rs) 



A 



(9) 



Equation (|9]l has a number of self-consistent solutions. For 
sub-ohmic baths one of these solutions is always A = 0. This 
solution corresponds to the complete localisation of the TLS 
[M ~ ±1) and is made self-consistent by the infra-red di- 
vergence of the integrand in Eq. (jSj for s < 1 and A = 0. 
This divergence arises from the divergence of the boson num- 
ber of the low frequency modes when subject to a static force, 
which causes \(f>±) to become orthogonal , 22 , 23 1 . As 
discussed in Refs. 1221 12311 . for sufficiently small a there are 
also finite solutions for A, and the physcially relevant one can 
be can be expressed analytically in terms of the Lambert W 
function l24ll . With this solution the groundstate energy can 
be written as a function of just the original system parameters 
and the magnetisation. To obtain the ground state magnetisa- 
tion we then simply minimise the ground state energy of Eq. 
^ w.rt. AI for fixed a, A and ujc- 
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GROUND STATE ENERGY, CRITICAL EXPONENTS AND 
CRITICAL COUPLINGS 

The ground state energy can be Taylor-expanded about 
M = 0. When this is done we find that for small M 
the energy takes the Ginzburg-Landau form E = co(a) + 
ci{a)M'^ + C2{a)M'^ + 0{M^) where Ci(a) are constants for 
fixed a,a;c, A. This form for the ground state energy guar- 
antees a second-order magnetic transition; above a critical 
coupling Q!c, a magnetisation grows continuously with M oc 
I a — ad 2 and the magnetic susceptibility x ^ \ct — ac\~^- 
The typical development of the continuous QPT is shown in 
Fig. [U The critical coupling ac is the coupling for which 
ci{oic) = 0. In the scaling limit this equation can be solved 
analytically and the critical coupling is obtained as, 



A, 



sin(7rs) 
27r(l - s) \ ~e 



(10) 



In order to obtain ac as a function of the original system pa- 
rameters, the renormalised tunneling matrix element at the 
critical point Ac must also be determined self-consistently 
from Eq. (|9]l- Fortunately, in the scaling limit Ac can also 
be found analytically, leading to the final prediction. 



sin(7rs)e"-'/2 f A 
27r(l - s) 



(11) 



Figure |2] shows the dependence of the critical couplings ob- 
tained from the variational ground state as a function of s for 
various values of utc and A = 1. Note that the predicted 
values agree well with those obtained by recent numerical 
QMC, NRG and SPSA studies, and reproduces the scaling 

ttc oc f ^' j previous seen in NRG and other approaches 
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FIG. 3. Magnetisation M as a function of {a — ac) /ac for a > ac, 
A = 1 and ujc = 10. For visibility, the curves have been multiplied 
by 1, 2, 4, 8 for s = 0.1, 0.2, 0.3, 0.4 respectively. Dashed lines cor- 
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FIG. 4. Expectation value (ax) as a function of a for A = 1 and 

UJc = 10. 



OBSERVABLES 
Magnetization and coherence 

Using the values of ac given in Eq. (fTTl i. we show in Fig. 
Q that the magnetisation obtained from the energy minimiza- 
tion does indeed behave like M oc (a — ac)^ close to the 
transition. The magnetisation data are a gain in good agree- 
ment with the QMC and SPSA results lH [ntl. Figure H 
shows the behaviour of the cohrence {a^) as a function of a 
for A/lUc = 0.1. We find that (Tx is always continuous at the 



transition, but 



iQ,^ is discontinuous. In the scaling limit 
the variational theory predicts {(Tx) = e~ ^(i-s) the critical 
point, and is thus independent of A/cjc- Above the transition 
{ax) decays faster, but persists well into the localised phase. 
The persistence and s-dependence of {<Tx) at and above the 
transition point seems to be generally consistent with the dy- 
namical NRG study of Anders et al. 1251], where it was found 



that a coherent tunneling peak in the equilibrium correlation 
function survives well into the localised phase for < s < 0.5 
and becomes more robust as s — >■ 0. This behaviour was also 
observed by Anders et al. in the non-equilibirum dynamics, 
where oscillation amplitudes at criticality become stronger as 
s descreases. This trend is described in our approach by the 
increasing value of (ax) around the transition as s decreases. 



Entanglement 

In Fig|5]the entanglement between the spin and the bosonic 
environment is plotted as a function of a. For globally pure 
states the entanglement £ is quantified by the entropy of 
the reduced density matrix of the TLS, and is thus given 
b y g = -p+ \o g{p+) - log(p_), where p± = i(l + 
V {o'x)'^ + (cz)^). We observe a very similar cusp-like maxi- 
mum in the entanglement £Max at ac to that found in a recent 
NRG study 1,261 1. As £Max is determined only by the value of 
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lution for A by prevent ing the infra-red divegence of the in- 
tegrand in Eq. ^ 11221 12311 . Despite the different frequency 
response of A and NA-modes, all modes are correlated in 
the same quahtative way with the TLS states = —fi-). 
Looking at Eq. (|2]i, one sees that with these dressing corre- 
lations, the ground state is not separable into TLS and bath 
variables, and thus there is entanglement between the TLS 
and bath at finite a. As dressing correlations monotonically 
suppress A as a is increased, the entanglement also increases 
monotonically as the transition is approached from the delo- 
calised phase. 



FIG. 5. Entanglement between the spin and bosonic environment as 
a function of a for A = 1 amd ojc = 10. 



Localised phase 



(ox) at OLc, we also predict that goes to a constant in the 
scaling limit with p± = 0.5(1 ± Ac/A). The entanglement 
quantifies the degree of non-classical correlations between the 
TLS and bath, and we shall now discuss the intutive physcial 
picture of the transition which arises from a microscopic anal- 
ysis of the correlations quantitfied by the entanglement. 



DISCUSSION 

Delocalised phase 

In the delocalised phase (A/ ~ 0) we find that /;+ = — /(- 
for all modes, and the AOD state thus coincides with the SH 
variational ansatz. The sub-Ohmic SH state was previously 
studied in Refs. 1221 12311 . and was shown to possess a tran- 
sition where A jumped discontinuously from a finite value to 
zero at a critical coupling cic- A similar discontinuous transi- 



tion was also found by flow-equation analysis ll27ll . However 
the discontinuous transition is in fact an artefact of the inflex- 
ibility of the SH ansatz, as it can only describe a delocalised 
phase or a completely localised (A = 0, M = 1) state; the 
discontinuous transition arises simply from the abrupt change 
from a delocalised to a completely localized state when the 
latter becomes energetically favourable. 

In the delocalised phase, the physics of the AOD/SH state is 
determined only by the renormalisation of A by the bath. The 
variational solution separates the bath into adiabatic modes 
(A-modes) and non-adiabatic modes (NA-modes) which have 
very diferrent frequency responses to the renormalised TLS 
tunneling. For the fast, high frequency A-modes (uji ^ A), 
the TLS tunneling appears to be a very slowly varying force, 
and the A-modes can adibatically adjust their displacements 
to maximise their interaction energy with the TLS (ji± ~ 
±5/W;"^) 11, El- The slow, NA-modes with cj/ < A can- 
not respond fast enough to follow the tunneling and their dis- 
placement is supressed at low frequency {ji± ~ ztgr^ajjA"^) 
I 22I I23II . From Eq ([8]l, one can see that this supressed dis- 
placement of NA-modes self-consistently permits a finite so- 



Above the transition (A/ 7^ 0), we find that a new energy 
scale appears in the problem, and the displacements /; j- are 
no longer of equal magnitude and opposite sign. Modes with 
uji ^ A/ A(l — A/^)^ 2 continue to be adiabatic, whilst non- 
adiabatic modes with ^ Af A(l — A/^)^2 now have dis- 
placements which have the same sign and grow at low fre- 
quency, //+ = //_ = -Mgiu:~^. 

Because the non-adiabatic mode displacements have the 
same sign, they are not correlated with the state of the TLS. 
As a result, the state of the system is essentially separable, and 
takes the form of a product state of NA-modes and the corre- 
lated TLS-A-mode wave function. The NA-mode wave func- 
tion has a finite displacement, and this appears to the TLS as 
an effective 'mean-field' -like magnetic field in the z-direction. 
Indeed, the displacement of the NA modes induced by the 
magnetization of the spin is precisely what one would obtain 
from a mean field treatment of the problem l [l^ . While the 
total average displacement of NA modes is finite, there is an 
infra-red divergence of the occupation number of these modes, 

as X); {A^i) NA oc M"^ dLuuj^~'^. However, from Eq. (|5]l, 
we see that the fact that fi± have the same sign leads to no su- 
pression of A from NA modes. Thus in spite of the infra-red 
divergence of the low frequency boson number, the coherence 
of the gound state remains finite at the transition and allows 
AI to grow continuously above the transition. 

As the magnetisation increases, the energy scale that di- 
vides the NA and A modes also increases, and more of the 
bath modes become non-adiabatic and uncorelated with the 
TLS. This trend is reflected in the monotonic descrease of 
entanglement above the transition, and also in the qualitative 
change in the oscillator-induced suppression of the TLS co- 
herence. In the localised phase, the suppression of {ax) be- 
comes increasingly due to the NA-mode bias, which caues the 
efective magnetic field seen by the TLS to point away from 
the X axis. In fact, as a 00 the suppression of A due to 
dressing vanishes, and the suppression of (ctx) is determined 
solely by the rotation of the ground state to lie along the ef- 
fective NA-mode magnetic field. 
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FIG. 6. Fidelity of the variational ansatz with the ground state de- 
termined by DMRG as a function of a for A = 1 amd ljc = 10. 
The DMRG simulation used 100 sites, Nt = 15 and 10 schmidt co- 
efficients were retained. At the critical coupling, the fidelity drops 
suddenly from finite values to zero. The inset shows the typically 
dramatic descrease in the fidelity with system size above the transi- 
tion. Inset data coiTesponds to s = 0.3, a — 0.3280 (ofc = 0.3267). 



FIDELITY OF AOD ANSATZ WITH DMRG GROUND 
STATE 



While there is general agreement between the observables 
obtained from the AOD ansatz and previous numerical studies 
of the transition, it is not posssible in most cases to compute 
the fidelity of AOD ansatz with the ground states obtained by 
these methods. In order to explore how close the AOD ansatz 
is to the true ground state, we simulated the ground state us- 
ing DMRG methods. When the results converge, the DMRG 
method produces an extremely accurate matrix product state 
representation of the ground state wavefunction of the system 
and bath, and the fidelity (overlap) of the AOD ansatz and 
the MPS ground state can be easily computed. To implement 
the DMRG simulation, the SBM is first mapped onto a semi- 
infinite chain of harmonic oscillators using an analytical uni- 
tary transformation previously used to study open-system dy- 
namics with time-adaptive DMRG (t-DMRG) (Si- Unlike the 
mapping onto a Wilson chain that is used in NRG approaches, 
our mapping does not require any discretisation of the spec- 
tral density of the bath, although our mapping can also be car- 
ried out analytically for linear and logarithmically-discretised 
baths, as shown in Ref. lfl9ll . Once mapped to a chain, the 
ground state of the model is obtained using imaginary-time 
t-DMRG evolution. Figure shows some representative re- 
sults for the fidelity F ~ (^'1$) between the AOD ansatz and 
the DMRG ground state | $) . In the delocalised phase we find 
that the DMRG results converge very rapidly, and that only 
the first few sites of the chain representation of the bath are 
appreciably excited. The fidelity with the AOD state (equiva- 
lent to the SH state) is extremely high considering the many- 
body nature of the system, and we conclude that the AOD/SH 
ansatz works well at weak coupling. As a approaches ac, the 
fidelity decreases considerably and is in fact zero above the 



transition. We believe that this suppression of F may arise 
from the truncation of bosonic fock space used in the DMRG 
simulation and the finite size of the harmonic chain. 

Applying the unitary transformation which maps the SBM 
to a harmonic chain to the AOD ansatz, one finds that that 
the average occupation of site n in the chain Nav{n) de- 
creases very rapidly in the delocalised phase, and an accu- 
rate simulation of the ground state with a small Ni, and chain 
length is possible. However as a approaches ac, the decay of 
Nav (n) along the chain becomes much weaker and contribu- 
tions from many sites need to be included. Although Nav {n) 
is bounded along the chain, the convergence w.r.t chain length 
in the DMRG simulations becomes much slower just below 
the transition, and for a fixed number of sites we observe a 
continuous cross-over rather than a sharp transition as a is in- 
creased. The magnetisation typically begins to appear in the 
DMRG simulation at a slightly below the ac predicted from 
the AOD ansatz, with the deviation decreasing as the number 
of chain sites is increased. This leads to the reduction of the 
fidelity seen in Fig. (|6]l as the transition is approached from 
below. Above the ac predicted from the AOD, the fidelity 
between DMRG and AOD groundstates becomes extremely 
small. This arises because Nav{n) actually diverges along the 
chain in the magnetic phase. This can be shown directly using 
the chain mapping 11911 . which predicts that Navin) oc n^"^'* 
as n — > oo when M is finite. As in NRG, the DMRG ap- 
proach uses a finite number of bosonic fock states Ni, to rep- 
resent each oscillator of the chain, and the divergence of the 
chain populations in the AOD ansatz for s < 0.5 cannot be 
described in this truncated basis. Consequently, the overlap 
(^'10) goes to zero exponentially with chain length, as the 
norm of the projection of the AOD ansatz onto the truncated 
fock space vanishes with increasing chain size. This rapid 
suppression of F with chain size is shown in the inset of Fig. 

Although of a different nature, the use of a finite Nb in the 
NRG description of the transition in < s < 0.5 also appears 
to cause problems, as discussed in Refs. lfl4[[l5ll . 

The result that Nav{n) oc n^^^'^ suggests that the cur- 
rent application of DMRG is unable to decribe the localised 
phase for s < 1/2. However, the AOD ansatz may provide 
a means to resolve this issue. The divergence of Nav (n) is 
ultimately related to the infra-red divergence, which is itself 
a consequence of the appearance of the 'mean field' displace- 
ment of the oscillators above the transition. The diverging 
occupations arising from this displacement could be removed 
by a simple unitary transformation H UHU^^, where 

U ~ cxp (^^M J2i 9i'^i^{o-i ^ ■ This corresponds to us- 
ing a displaced oscillator basis where the effects of the mean- 
field displacement (computed from the AOD ansatz) are re- 
moved from the problem. Performing this unitary tranforam- 
tion and then mapping the resulting SBM onto a chain, we 
find that Nav{n) computed from the AOD ansatz now de- 
creases along the chain at a similar rate to that found in the 
delocalised phase. Therefore, using the AOD prediction for 
M and a unitary transormation may allow an efficient sim- 
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ulation of the localised phase using DMRG in a truncated 
bosonic fock space. In the likely event that the AOD mag- 
netisation is not exactly the true ground state magnetisation 
of the model, the unitary transformation above will not com- 
pletely cancel out the infra-red divergence and the siumlation 
will converge slowly. Therefore, by varying the value of M 
used in the unitary transformation, we can seek out the true 
value by looking for the value with the fastest convergence 
time. This idea of removing the effect of the mean-field dis- 
placment was orginally used in the SPSA study where it was 
found to yield an efficient and accurate simulation of the mag- 
netic phase UtIi . However, this method will only be applica- 
ble to the case described here (s < 0.5), where a mean-field 
bias arises from the bath. For s > 0.5 the transition is not 
expected to be mean-field, and it is perhaps interesting that in 
our chain mapping the asymptotic values of Nav (n) are qual- 
itatively different above and below s = 0.5. 



CONCLUSION 

We have presented a strongly-correlated variational ansatz 
for the ground state of the sub-Ohmic spin-boson model and 
have found that it undergoes a continuous quantum phase tran- 
sition at a critical coupling strength to the bath. The critical 
properties and the values of observables that it predicts are 
very good approximations to those found by numerical meth- 
ods, and by an analysis of the microscopic structure of the 
ansatz we have been able to describe the rich physics of the 
adiabatic and non-adiabatic modes that may drive this transi- 
tion. We have also found that an infra-red divegence of the low 
frequency bosons always accompanies the transition in this 
theory, and this may be of relevance for numerical approaches 
to this problem, as illustrated with our numerical DMRG sim- 
ulation of the sub-Ohmic ground state. 
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